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Abstract 



In this paper, using the Newman-Penrose formahsm, we find the Maxweh 



equations in NUT space and after separation into angular and radial compo- 



nents solve them analytically. All the angular equations are solved in terms 



of Jaccobi polynomials. The radial equations are transformed into Hyperge- 




ometric and Heun's equations with the right hand sides including terms of 



the expansion of this parameter. 



different order in the frequency of the perturbation which allow solutions in 



PACS No. 04.20.-q 



Typeset using REVTeX 



* Email: nouri@khayam.ut.ac.ir 



1 



I. INTRODUCTION 



An effective way to understand and analyse characteristics of an spacetime is to study 
its behaviour under different kind of perturbations including the most studied one electro- 
magnetic perturbations. Indeed interest in the analytic solutions of Maxwell's equations in 
curved spacetimes, specially in the well known solutions of Einstein field equations, arose 
from the study of stability of those spacetimes under electromagnetic perturbations. On 
the other hand the (Taub)-NUT solution of Einstein vacuum field equations [1,2] has found 
importance both in the discussions on magnetic monopoles in gauge theories [3-4] and in 
string theory [5-6], though by itself interpreted as the metric of a mass M endowed with a 
gravitomagnetic monopole with strength /, the so called NUT parameter or magnetic mass 
[7]. Therefore it seems appropriate to study electromagnetic waves in NUT space through 
the solutions of the Maxwell equations in this spacetime and use the results to obtain a 
deeper understanding of the physical aspects of NUT space . 

In a recent paper [8] perturbations of the Kerr-Taub-NUT spacetime by massless fields of 
spin s < 2 are studied by introducing a master equation describing these fields. The au- 
thors have also discussed the coupling between the gravitomagnetic monopole moment / and 
the spin of the perturbation field in the context of gravitomagnetism. In [9] Klein-Gordon 
and Dirac equations in Kerr-NUT space have been studied with an emphasis on a duality 
transformation involving the exchange of mass and NUT parameters. It is shown that the 
same duality transformation also exchanges the scalar and spinor perturbations. In the 
present article, as stated in the abstract, we have a more restricted task in front of us and 
that is to study electromagnetic perturbations (massless spin 1 field) in NUT space (and not 
Taub-NUT space) but we intend here to solve the resulted equations more or less anlytically. 

To do so we start with the fact that as in the case of the Kerr metric, because NUT 
space is stationary and axisymmetric ^, one could express a general solution to the Maxwell 



For a review on the physical symmetries of NUT space refer to [7] and [10]. 
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equations in this space as a superposition of different modes with a time- and (^dependence 
given by [11] 

where m = 0,±l,±2,.... and to is the frequency of each mode. Treating the electromagnetic 
field as a perturbation, it is natural to think of c<; as a small parameter in terms of which 
the appropriate functions could be expanded if necessary. 



II. MAXWELL EQUATIONS IN NUT SPACE 

In this section using Newman- Penrose formalism [12] we separate the remaining two 
variables, r and 9 and obtain the radial and angular equations involved. To do so we start 
with the NUT metric given by the line element 

ds^ = f{r) {dt - 2lcos9d<pf - -T^dr^ - (r^ + l^){d9^ + sin^9df) (2) 

where /(r) = 1 — ^ > and I is called the magnetic mass or NUT factor. Using the 

tangent vectors to the null geodesies in NUT space, one can easily show that the following 
set of null vectors form a null-tetrad basis for the above metric 

l^^{^ylM , I,^{l,-^y0,-2lcos9) 
= 0, 0) , n,= (^, 1 0, -/cos^/(r)) 



1 7-2 + P 

= ^= (2i/cotg,0,l,icscg) , = \ -(0, 0, -1, -isin^) 

^2(r2 + P) V 2 



1 /r^ + P 

fh"^ (-2i/cotg,0, l,-icsc6>) , % = V — ^ — (0, 0, -l,isin^) (3) 

/2(r2 + /2) V 2 



Using the above null tetrad one can find the following spin coefficients of NUT metric; 

«; = 0, (7 = 0, A = 0, z^ = 0, r = 0, 7r = 

_—r^-\l _ f{r) -r + il 

^~ r^ + P ' ^" ^V2 + /2^ 
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■ 7 = 7 / (^) + 



2(r2 + P) ' ' 4 V r2 + 

—1 cot^ ^ 1 cot^ 

« = -7^ — , ; d = - , = (4) 

2 ^2(r2 + /2) 2^2(r2 + Z2) 

Having found the spin coefficients of NUT metric we can now write the Maxwell equations in 

this space using the Newman-Penrose formalism [12]. For this we replace the antisymmetric 

Maxwell tensor F^,^ with the following three compiea; scalars [11]; 

02 = F42 = F^.m'^n'^ (5) 

where Fij — 1,2,3,4) and F^i^ (//, z/ = 0,1,2,3) are the components of the Maxwell 
tensor in the tetrad and tensor bases respectively. One can invert equations (5) to find 

F^u = 2 4)omii^n^] + 0i(n[^/^] + mi^rh^]) + 02^[M"^i/] +C . C 

In terms of the above scalars the Maxwell equations become; 

D(j)i - 5*00 = (tt - 2q;)0o + 2p0i - fi;02 

D02 - 5*01 = -A0O + 27r0i + (p - 2e)02 
501 - P0O = (a* - 27)00 + 2r0i - (702 

502 - 7^01 = -z/00 + 2/X01 + (r - 2/3)02 (6) 

Where in the above equations D,V,S and 5* are the special symbols for the basis vectors 
1, n, m and m when they are considered as the directional derivatives. For NUT space these 
are; 

~ ~ /(r) dt ^ dr 

^ 15 f d 
n = V = -— - - — 
2dt 2dr 
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1 /^., ^9 d id 
^ = S^ , (2i/cot^- + — + 



^J2(r^~+P) 9t do smO 84) 

^^^*= I ^ (-2i/cotg|- + ^ - ^|-) (7) 
/2(r2 + p) dt 89 sin9d(f) 

Now if wc use the general form (1) for the {t, <p)- dependence of the wave, the above operators 

take the following forms; 



2(r2 + P) 



2 2"^ 2 

—2lujcoW -\- de — mcsc9) 



S* = , ^ (2lujcot0 + dg + mcsc9) (8) 
^2(r2 + /2)^ 

Now we define the above operators in the following way; 

m. — S= , ^ =>^o , m — S*— , ^ =^0 

^2(r2 + /2) y'2(r2 + /2) 

where 

= Or + iTTT + ^'^ — X — 
/(r) A 

w o . ^ r-M 

'^n = (^e> + Q + ncot9 

Cl = de-Q + ncotO (9) 

in which A = — 2Mr — P and Q = 2lujcot9 + mcsc9. We note that is (up to a 
factor / \ ^ ) nothing but one of the weighted derivative operators introduced in the 

y 2(r2+/2) 

compacted-spin coefficient formalism [13]. Now substituting the above operators in the 
Maxwell equations (6) we have; 

^2(r2 + P) + 



1 r 

£o01 = {1^0 + n , ,o )02 



1 



^2(r2 + P) 
1 



401 = 



fir) 



-{Vl 



'2{r^ + P) 2 r^ + P' 

The above equations will look more compact if we change the variables to; 



where now, using the definition of the operator we have; 



r'^ + l 



r^ + P 



4$2 = -A{vi + 



4$i = -A{V\ 



r^ + P' 
r 

r^ + P 



(10) 

(11) 
(12) 
(13) 



The commutativity of the operators Cq and — A('Do + ^sq^) enables us to ehminate $i 
from equations (10) and (13) and get; 



r^ + P' 



(14) 



where we have used the fact that VnA — AVn+i- Similarly the commutativity of the 
operators £0 and —A{Vq + ^q^) will allow us to eliminate $1 from equations (11) and (12) 
and get; 

$2(r,^) = (15) 



£o4 + A(pS + ^^)(^^o 



r^ + P'' " r^ + P' 

Now our task is to find solutions to the equations (14) and (15) by separating them into 
angular and radial parts. This separation can be achieved by choosing the following general 
form for ^Q{r, 9) and $2(^", ^); 



$o(r, 9) = no{r)Qo{9) and $2(r, 9) = 7^2(r)e2(^) 



(16) 



Substituting (16) in (14) and (15) we have; 

4£ieo(^) = {de -Q){d0 + Q + cote)eo{0) = -COoiO) (17a) 

^(^1 + - :;:^2)^o{r) = e7^o (176) 

£o402(^?) = {de + g)(ae - Q + cot^)e2(^^) = -^e2{0) (18a) 

where ^ is the separation constant. These are analogous to the Teukolsky equations in Kerr 
space [11,14]. Note that we have not distinguished between the separation constants that 
derived from equations (14) and (15) and the reason is that equations (17a) and (18a) deter- 
mine the same set of eigenvalues ^. This can be seen by considering equation (17a) where, 
apart from the regularity requirement of ©o(^) at ^ = and ^ = tt for the determination 
of ^, the operator acting on 02 (^) in equation (18a) is the same as the operator acting 
on Go in (17a) if we replace 9 hy n — 9 and this is the case because the general relation 
^n{9) — —Cl{7r — 9) holds. So a proper solution , ©o(^;0) of equation (17a) with eigenvalue 
^ , is also a solution of equation (18a) with the same eigenvalue, if we replace ^ by tt — ^ in 

Qo{9;0- 



III. SOLUTION TO THE ANGULAR EQUATION 

In this section we solve equation (17a) analytically. It is obvious from the form of the 
equation (18a) and the discussion above that it will have the same set of eigenvalues as 
(17a). First of all we rewrite equation (17a) in the following form; 

{dl + cot9d0)eo{9) + eo{9)d0{Q + cot9) - {Q^ + Qcot9 - 0^o{9) = 

Substituting for Q from its definition, the above equation becomes; 

{dl + cot9de)eo{9) - -^{21lv + l)eo{9)- 

sm u 

J— + 4:lWcot^9 + — J- + Alum — j- + 2lujcot^9 - ^ QoiO) = 

y sin y sin sin 6' J 
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Now changing the variable to x — cos9 we have; 



(1-^)372-2^:1- 



d {2mx + Ab'^x^ + + Ahmx + 2hx'^ + 26+1) 



dx 



dx 



Qo{x) = (19) 



where b — lu. Dominant behaviour of the eigenfunctions of the above equation at singular 
points a; = ±1 are given by (1 + x)^~^ respectively. Introducing a — |1 + 26 + m| and 

/3 = |1 + 26 — m|, and substituting for Bq the following expression; 

i^-^r/\'-±^f-u{x) 



in (19) it will be transformed into; 

(l-x^)^^ + [(3-a-{2 + a + (3)x 
dx 



dU{x 



+[e-(^)(^ + l)+46^ + 26][/(a;) = 

(20a) 



This is the differential equation satisfied by the Jacobi polynomials [15,16] 



:-ir 



2"n! 



dx^' 



if we identify; 



^ + 46^ + 26=[n+(^)][n+(^ + l)]=j(j + l) 
where j = n + {^-k^)- So 0o could be written explicitly as follows 



Qo{x) = J , (1 - x)-t(l + x)-^^ [{I - a;)"+"(l + ^r^-" 
2" 2 2 77,! cia: 



d" 



(206) 



(20c) 



IV. SOLUTIONS TO THE RADIAL EQUATIONS 

Now we try to solve the radial equations (17b) and (18b) by replacing for the operators 
V and Vi and their complex conjugates from (9); 



^2 _|_ 72 ' 



A2 



/ 2-2ia;(r-iO-g 4i/(r - M) 2ilr - P \ ^ 
[ A A(r2 + Z2) + (^2 + ^2)2 ) 



(21a) 
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_l_ 



+ + A + (r2 + J ^2 - (216) 

where A — r'^ — 2Mr — P — (r — r+)(r — r_) with r± = M ± V + P and r+ being the 
position of the horizon. Therefore in our discussion below we always have r > r+. 

A. Solution to the equation (21b) 

We start by solving equation (21b) which looks simpler. This equation has four regular 
singularities at r = r± and r = ±i/ and an irregular singularity at r = oo. To remove the 
singularities at r = ±iZ we use the following transformation ; 

^2(r) = (^)'/'V2(r) (22a) 
so that equation (21b) becomes (after some manipulation); 



where 



(226) 

A = r+(2Mu^ - iou) + 6^ 
S = r_(2Ma;2-iu;) + 62 
C = 2M{2Muj^ + iuj) + 2(6^ + 6) - ^ 

D = 2Mu'^ + iu (23) 



Now if we define 



r — r+ = —ex & r — r_ = £(1 — x) (24) 
where e = 2(M2 + ^2)1/2 and take 

V2 = (-x)"(l - xfV2 
9 



(where a and P are constants to be determined later) then equation (22b) will transform 
into; 



A 

— + 



B 



C 



'1-x) 



x(l — x) 



V2ix) 



+ 



eD( 



1 



r „ „ a(a-l) 2a(3 (3(l3 - 1) 
'+e^u}^ + ^^ — 7-—^ + 



V2ix) = 



(25) 



1 — xx'' ' x^ x{l — x) {1 — x)'^ 

Now in order to eliminate the terms proportional to ^ and (jz^, ck and (3 should take one 



of the following values. ; 



1 



/3± = -(i±vr^^) 



(26) 



As in the case of the Kerr spacetime one can show that different choices of the above param- 
eters correspond to solutions which satisfy the incoming and outgoing boundary conditions 
[17]. For example choosing the pair (q!_,/3_|_) corresponds to the solution satisfying the 
incoming boundary condition^. With the above choice of a and /3 equation (25) reduces to; 

1 



dl%{x) + 2{^--^)d^V,{x) + 



eD( 



1 — X X 



c 



I 2 2 

+ e to 



2ap 



Voix) = 



(27) 



x(l — x) x{l — x) 

Now choosing 1^2(3:^) = exp{ieux)Vi{x) the above equation can be written in the following 
form; 

x{l - x)dlV2{x) + [2a - 2{a + P)x]d^V2{x) - pqV2{x) = 

{{2M + e){2Mu'^ + iu) + 2{3lW + 2luj) - j{j + 1) + 2aP - 2ieua)V2{x) + 

{2isij[a + /9] - 2e{2Muj^ + iuj))xV2{x) - 2ieujx{l - x)d^V2{x) - pqV2{x) (28a) 

where we have substituted for C, D and ^ from equations (20a) and (23) and added —pqV2{x) 
on both sides to make the left hand side of the equation to look like the hypergeometric 



^This could also be seen by the fact that a = limit of the Kerr case and / = limit of NUT case 
should both reduce to the electromagnetic perturbation of schwarzschild spacetime. 
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equation with the identification 

p + q^2(a + /3)-l (28b) 

The right hand side of the above equation includes terms of different order in to and therefore 
it is suitable for finding the solution in the expansion of this parameter and it is obvious 
that the zero-th order solutions of the equation are the hypergeometric series provided we 
choose 

pq^2aP-j{j + l) (28c) 

To write the zero-th order solution more explicitly we note that for 2a a non-integer the 
general solution in terms of hypergeometric series could be written as follows; 

V2(x) = AF{p, q- 2a- x) + Bx^~'^'^F{p - 2a + l,q - 2a + 1;2 - 2a; x) 

hence to zeroth order in to we have; 

r — ilir — r+ r — r^g 



7^2 r 



r + il e e 
(aF{p, q; 2a; !±^) + B{^^^^Y-^''F{p - 2a + l,q - 2a + 1;2 - 2a; (28c/) 



where A and B are constants and p and q are solutions of equations (28b) and (28c). 

B. Solution to the equation (21a) 

Amazingly enough one can see that the same procedure used to transform equation (21b) 
to (28) can also be applied to equation (21a). The reason for this is the fact that the same 
transformation as in (22a) will remove the singularities of (21a) at r = ±i/, transforming it 
into; 

r — r+ r — 

/A' R' C 1 1 \ 

^ D'{^— + ^—)+u;']Voir)^0 (29) 



(r — r+)2 (r — r_)2 {r — r+){r — r-) r — r+ r — r-' 
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where 



^o(r) = C—^y^'Mr) (30) 



and 

B' = r_(2Ma;2 + ia;) + 6^ 

C = 2M{2Mu'^ - iu) + 2(6^ + 6) - ^' 

D' = 2Ma;^ - io; (31) 

It is notable that these equations are the same as in (23) with iuj transformed into —iui and 
^ to ^' = ^ — 2. Now applying the same change of the variable as in (24) and substituting 
for 

Vo(r) = (-x)"'+^(l - xf^^ ex^{ieujx)Vo{x) (32o) 

where 

a'± = i(-3± V1-4A') , /3± = ^(-3± Vl-45') (326) 
equation (29) will transform into 

x{l - x)dlVo{x) + [2a' - 2{a' + ^')x]dJi?o{x) - p'q'Vo{x) = 
[2i£io{a' + /?') + 2e(2M£x;^ - ia;)]xVo(x) + 2iec<;x(l - x)d^Vo [x) — p'q'Vo{x) + 
(2M + £)(2Ma;^ - iu) + 2{3lW + 2luj) - + 1) + 2{a' + l){p' + 1) + 2i£a;Q;'] Vo{x) (33) 

where 

/(/ + l) = e' + 46' + 26 (34) 

As one can see the left hand side of the above equation is also the hypergeometric equation 
with the identification p' + q' — 2{a' + /?') — 1. As in equation (28a) the right hand side of 
this equation includes terms of different order in u> and it is therefore suitable for finding 
the solution in the expansion of this parameter. Again the zero-th order solution is the 
hypergeometric function provided we choose 

pV = /(/ + 1)-2(«' + 1)(/3' + 1). 
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V. EQUATION SATISFIED BY $i AND ITS SOLUTION 



To complete our discussion we need to solve the equation satisfied by scalar 0i which, 
using equations (11) and (12), can be written in the following form; 



$i(r,^) = (35) 



4^o + A(I),--^)(pt + r_^) 
Now upon choosing the form 

$l(r,^) = 7^l(r)el 
will transform into the following angular and radial equations; 

4>Co©i(^) ^{de-Q + cote){de + Q)e,ie) = -rjOiie) (36) 



^(^1 - ;^)(^o + = (37) 

Now our task is to solve the above two equations. 



A. solution to the angular equation 

We start with (36) which can be solved in the same way as we solved equations (17a) 
and (18a). Substituting for Q from its definition and changing the variable to a; = cos^, 
equation (36) will transform into; 



2, n ^ d (2mx + Ah^x^ + w? + Ahmx + 2hx^ + 26) 

1 - ^ - 2a^3- - ^ 2 + ^ 

dx dx 1 — 



61(0;) = (38) 



where h = Ico. One can see that the only difference between this equation and equation (19) 
is the disappearance of 1 in the numerator of 1 — x^. In a similar way to equation (19), the 
dominant behavior of the eigenfunctions of the above equation at singular points a; = ±1 
are given by (1 T^;) ^ where B± = |1 + 26±m|. Introducing a' = {B+ + 1)^/^(,B+ — 1)^/^ 
and P' = {B- + lY^'^{B- — 1)^^^, and substituting for 0i the following expression; 
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in (38) it will be transformed into; 

(39) 

This is again the differential equation satisfied by the Jacobi polynomials if we identify; 

77 + 46^ + 26 = [n + (^^)][n + (^^ + 1)] = j{j + 1) (40) 
where j = n + )■ 

B. solution to the radial equation 

Substituting for the operators from their definitions given in (9), the radial equation (37) 
can be written in the following form; 

^UAr) + 2(^ - ^)OrnAr) + {f^^) n^r) 

One can remove the singularity a.t r = il by the following transformation; 
so that equation (41) transforms into; 

(43, 

Then by the following change of variable, 

r — r+ 
2M 

the above equation transforms into 



= —z 



1 1 2 \ apz-q [(^r+ - 2Mzf + P^^ 



^ z-1 z-aj z{z-l){z-a) 2M^ z^{z-lf ^ ' 
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where 



and 



a + (3^-1 ; a/3=^{2ul + r]) (446) 



q — a 3a : a = — — (44c) 

Comparing equations (44a-c) with Heun's differential equation [17] defined by 

dz^ \z z — 1 z — a J dz z{z — l){z — a) 

where 

a ^ 0,1 and 7 + 5 + e = Q; + /3 + l 

we find out that equation (44a), along with the fact that \a\ > 1, forms the Heun equation 
with the right hand side including a term of second order in cu, allowing solutions in the 
expansion of this parameter. Again the zeroth order solutions are solutions to the Heun 
equation but before looking for that the following two points worth to be mentioned; 
First we note that despite the fact that q = aa/3, the left hand side of equation (44a) does 
not degenerate into hypergeometric equation because e ^ and secondly, due to the fact 
that all the parameters 7, 5 and e are integers we expect that the general solution to the 
Heun equation involve logarithmic terms. 

To write the explicit zero-th order solution to the equation (44a) we only consider the basic 
power series solution ^ to the Heun equation denoted by [18]; 

00 

Hl(a,q;a, (3,^,5; z) ^^CjZ^ , (co = 1) (46) 

j=0 

where HI stands for 'Heun local' solution and apart from the normalized coefficint Cq the 
other coefRcints are given by the following three-term recursion relation; 

—qco + 0701 = (47) 



^For detailed dicussion on possible different solutions to Heun's equation see [18-19] 
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Pjcj-i - {Qj + q)cj + Rjcj+i = , (j > 1) (48) 

in which ; 

Pj^{j-l + a){j-2-a) (49a) 
Qj=j[j{l + a) + a-2] (496) 
Rj = a{j + 1)2 (49c) 
where we used relations (44b) and the fact that in our case, left hand side of equation (44a), 

^ — 5 — 1 and e — —2 

Therefore the zero-th order solution to the equation (44a) is given by; 

It is also interesting to note that Heun's equation and its biconfluent extension also 
appear in the discussions on the electromagnetic perturbations of Kerr-de Sitter and Kasner 
spacetimes respectively and indeed furnish exact solutions to the Maxwell equations in these 
spacetimes [20,21]. 

VI. RESULTS AND DISCUSSION 

In this paper we have considered the Maxwell equations in NUT space and solved them 
analytically using the Newman- Penrose null tetrad formalism, which is best adapted for 
treating type-D spacetimes. It is shown that after separation of the equations, solutions to 
the angular equations could be given in terms of the Jacobi polynomials. The radial equa- 
tions on the other hand are transformed into hypergeometric and Heun's equations which on 
the right hand side have terms of different order in uj. This fact enables one to find solutions 
in the expansion of this parameter with the zero-th order solutions being the hypergeometric 
function and possible solutions to the Heun equation respectively. To compare the results 
with the Schwarzscild case we note that the appearance of the perturbation frequency uj in 
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the angular equations (19) and (38) , only through the parameter b — Ico, means that unlike 
Schwarzschild, in the stabihty discussions of NUT space, we can not only consider radial 
equations and angular equations should also be considered. In this respect the behaviour 
of NUT space under electromagnetic pertubations is very similar to that of Kerr in which 
angular equations depend on the frequency through the combination au [11] . So in both 
cases going over to the Schwarzschild space will remove the perturbation parameter from 
the angular equation and the stability discussion will only need a study of radial equations. 
Another point deserving attention is to check whether, as in the case of the Kerr geometry, 
it is possible to write the solutions to the radial equations in the form of an expansion in 
terms of the coulomb wave functions [17,22] . 
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